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ABSTRACT 

Q>^ ! A nonlocal and nonlinear theory of hadrons, equivalent to the color singlet sector 

rH ! two dimensional QCD, is constructed. The phase space space of this theory is an infinite 

dimensional Grassmannian. The baryon number of QCD corresponds to a topological 

^ ! invariant ('virtual rank') of the Grassmannian. It is shown that the hadron theory has 

^ . topological solitons corresponding to the baryons of QCD. -^ plays the role of h in this 

in 

theory; A''^ must be an integer for topological reasons. We also describe the quantization of 
--^ , a toy model with a finite dimensional Grassmannian as the phase space. In an appendix, 

Q\ ! we show that the usual Hartree-Fock theory of atomic and condensed matter physics has 

-)— > , a natural formulation in terms of infinite dimensional Grassmannians. 

Oh: 

> 
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1. Introduction 

Qunatum Chromodynamics (QCD) is the generally accepted theory of strong interac- 
tions. It describes the strongly interacting particles (hadrons) in terms of their constituents, 
the quarks and gluons. However the quarks and gluons themselves are not directly ob- 
servable; only their color singlet bound states (hadrons) exist as isolated particles. A 
challenging problem of particle physics is to construct a theory of strong interactions di- 
rectly in terms of hadrons. There are indications that such a theory exists and is a sort of 
string theory [1], although a complete understanding is not available yet. 

In this paper we will present such a theory of hadrons (Quantum Hadrondynamics, 
QHD) in two spatial space-time dimensions. It will be shown that the color singlet sector 
of two dimensional QCD (2DQCD) is equivalent to QHD for all energy scales and number 
of colors. At low energies the theory tends to a nonlinear sigma model. For large N^. it 
tends to a clasical theory; the quantum fluctuations in QHD are of order -^. This classical 
hadron theory however has no relation to classical chromodynamics; it is equivalent instead 
to the large Nc limit of quantum chromodynamics. 

QHD can be viewed as an extension of the ideas of Berezin [2], [3] to QCD. There 
is also some previous work by Kikkawa [4] in the same spirit, although our theory 
differs in some important aspects. One main point clarified by our work is the origin of 
baryon number as a homotopic invariant in QHD. For this, as well as other reasons, an 
understanding of the quadratic constraints satisfied by the field variable is crucial. The 
mathematical framework necessary is the theory of infinite dimensional Grassmannians as 
formulated in the book by Pressley and Segal [5] . A less explicit summary of this work 
was given in [6]. 

The main distinguishing feature of QHD is that it is not a local quantum field theory. 
The field variable depends on a pair of space- time points separated by a null distance. 



Also, the theory is highly nonlinear, the phase space of the theory being a curved manifold 
(the infinite dimensional Grassmannian). We will first describe the classical limit of QHD, 
which is equivalent to large Nc QCD. Small oscillations around the vacuum describe mesons 
and can be decribed by a linear approximation to large Nc QCD constructed in the work 
of 't Hooft [7]. Our approach will construct the complete theory, which has interactions 
at every possible order, since the field variable takes values in a curved manifold. It 
is important to note that the infinite Nc limit of QCD is equivalent to an interacting, 
highly nonlinear, theory. It is only if we further restrict ouselves to the small oscillations 
around the vacuum that it becomes a theory of free mesons. If large Nc QCD were a 
free field theory as is occassionally claimed there would be no soliton solutions, which are 
necessary to describe baryons as argued by Witten [8]. That baryons should be solitons 
of a theory whose small oscillations are mesons was proposed independently by Skyrme 
[9] before QCD. It is now known that this idea is consistent with QCD [10]. QHD has 
a homotopically conserved quantum number which can be identified with baryon number. 
There are static solutions to our theory carrying this quantum number, which do look like 
baryons [11]. 

Our classical hadron theory is equivalent to the sum over planar diagrams of 2DQCD. 
The field variable M{x, y) of our theory is the 'master field'; we will obtain its equations in 
closed form. It would have been too hard to obtain by actually summing diagrams, since 
there are vortices at every order. It is only the geometrical understanding of the phase 
space as a Grassmannian that makes the formulation of the theory possible. 

Quantization of our hadron theory corresponds to studying QCD at finite Nc- In the 
hadron theory, Nc plays the role of ^ ; it is required to be a positive integer for homotopic 
reasons. The solution of the hadron theory in the semiclassical approximation corresponds 
to the -^ expansion of 2DQCD. 



2. Classical Dynamics on Finite Dimensional Grassmannians 

2. 1. Parametrizations of the Grassmannian 

We are eventually interested in the infinite dimensional Grassmanian, which is the 
phase space of classical hadron dynamics. However, the particular infinite dimensional 
Grassmannian ( the 'restricted' Grassmannian of G. Segal [5]) we need, is very similar 
to the finite dimensional counterparts. It is useful therefore to study the analogous finite 
dimensional spaces first. This can be viewed as a 'regularization' of our theory. 

For a positive integer M we will define the Grassmannian Gtm to be the set of all 
M X M hermitean matrices satisfying a quadratic constraint: 

G'rM = {$|$^ = $;$^ = l}. (1) 

The eigenvalues of $ will then be ±1. Each point $ G Gtm picks out a subspace of C^, 
the eigenspace of $ with eigenvalue —1. Using this correspondence the Grassmannian can 
be viewed as the set of subspaces of C^ . This is the conventional definition [12]. 

The trace tr$ is an integer: 

tr$ = M - 2m (2) 

where m is the number of eigenvalues equal to —1. Thus tr$ is invariant under a continuous 
deformation of $. In fact Gtm is a union of components labelled bym = l,---M — 1. The 
case m = {m = M) is trivial, containing just one point $ = 1($ = —1). 

M 
Gtm = IJ Grra,M- (3) 

771 = 

Each component can be viewed as a coset space of the unitary group: 

Grm,M = U{M)/U{m)xU{M -m). (4) 
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To see this, note that any hermitean matrix can be diagonahzed by a unitary transforma- 
tion. There wiU be precisely m eigenvalues equal to —1 and M — m equal to +1. Thus, 
for each $ G Grm,M, there is a, g E U{M) such that 

$ = geg^ (5) 

where 

e = n 1 • ^^) 

\ U i-{M-m)x{M-m) J 

Both g and gh correspond to the same $ if /i commutes with e. The subgroup of elements 
U{M) that commutes with e is U{m) x U{M — m), consisting of unitary matrices that are 
block diagonal: 

-(';? :)■ (^) 

Thus there is a one-one correspondence between $ G Grm,M and the coset space 

U{M)/U{m) xU{M -m). 

We see now that each component Grm,M is a connected, compact manifold of dimen- 
sion 2m(M — m). Also the map $ ^ — $ is a diffeomorphism of Grm,M to GrM-m,M- 
If M = 2m, this is a map of Grm,M to itself. The group U{M) acts transitively on each 
connected component Grm,M by the action 

$ ^ (/^(/t. (8) 

The simplest special case of a Grassmannian is 

Gn,2 = U{2)/U{1) X U{l)SU{2)/U{l) = S\ (9) 

In fact the description of the general Grassmannian as the set of matrices satisfying a 
quadratic considtion is similiar to that of S'^ as the set of vectors of unit length. It is useful 
to keep this simple example in mind as we develop the general theory. More generally. 



Gfi,M is the complex projective space CP^ ^ each point of which is a one dimensional 
of subspace in C^ . 

A vector field on the Grassmannian en be thought of as a matrix valued function V^($) 
satisfying 

[y($),$]+ :=y($)$ + $y($) = 0. (10) 

(If $ = e, this means that V is off-diagonal ^ = I f n,)-) ^^y such matrix valued 
function is of the form 

y($) = [$,t/($)] (11) 

for some other function t/($). However, U is not uniquely determined by V: the trans- 
formation t/($) -^ t/($) + [$, A($)]+ wiU leave V^($) unchanged. Thus t/($) is a sort of 
'potential' for V^($), since the commutator with e is like a differentiation. (It satisfies the 
formal rules for cyclic cohomology; see [13]). If V^($) is a vector field, we have also the 
identity 

[$,[$, y($)] =4y($). (12) 

2. 2. Coordinates on the Grassmannian 

It is sometimes convenient to solve the constraint on $ by introducing explicit co- 
ordinates. Just as for 5"^, one needs several coordinate charts to cover all of the Grass- 
mannian. 

Consider first the case of Gr 1^2 = CP^ . The one dimensional subspace picked out by 
$ consists of vectors of type AZ\ i = 1,2. The vector Z"^ can be schosen to have unit 
length, so that it forms an 'orthonormal basis' in this one dimensional subspace. Since Z^ 
is an eigenvector of $ with eigenvalue —1 and the other eigenvalue is 1, 

^=5l-2ZtZ^. (13) 



Note that Z* -^ Z^h will leave $ invariant if h E U{1). The point $ = e corresponds to 
Z^ — I \ . More generally, we can choose the unit vector Z^ in the form 

' "' (14) 



v/(i + </>*</)) \4>* 

4> being some complex number. The magnitude of the first entry is fixed by the length of 
Z^ and its phase can be chosen as above by an appropriate choice of h. This gives, 

One can check easily that $^ = 1 and tr$ = 0. This coordinate system overs an open 
neighborhood of e, but it breaks down at $ = I j . Another co-ordinate system. 



based on the choice 



and 



^^^7(r^('' ^''^ 



is well-defined in a neighborhood of this point. The transformation that links these systems 



IS 

,, 1 



(18) 



Since the space can be covered by coordinate charts related by holomorphic(analytic) 
cordinate transformations, we see that Gr 1^2 is a one dimensional complex manifold. 

The above construction can be generalized to an arbitrary Grassmannian Grm,M- We 
will write 

^1=51- 2ZIZ^\ (19) 

where Z^"- for a = 1, ■ • ■ m is an orthonormal basis for the eigenspace of $ with eigenvalues 
—1. Z can be regarded as an M x to matrix and the orthonormality relation is 

ztz = 1 (20) 
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In matrix notation, $ = \ — 2ZZ^. Z and Zh give the same $ if /i G U{m). We can choose 
Z to be 

Z=( !\)[1 + <P<P'^]-'^' (21) 



where (p is an m x (M — m) complex matrix. (The 1 in the above equation stands for an 
m X m identity matrix. 1 + (jxj)'^ is a positive invertible m x m matrix so it has a weU- 
defined —1/2 power.) This gives 



+ 



io l) V(/'"^(l + </'</>■'■ )"^ ^t(i I ^^t^-l^ ] • ^ > 



One can check exphcitly that $^ = 1, tr$ = M — 2m. 

The action of U{M) on $, $ ^ 9^9^ can be seen to be the transformation 

(j)^go(j)={c + d(j)){a + h(i))~^ (23) 

, { a h^ 

where g = , 

Again this coordinate system breaks down at the points $ = Treyr"^, tt being a per- 
mutation that interchanges a negative with a positive eigenvalue of e. ( Permutation of a 
pir of eigenvectors of e is a unitary transformation in C^). If we define, 

(/), = TT O (/, (24) 

the new coordinate will remain well-defined even at $ = Treyr"^. Thus we can cover all 
of Grm,M by I J coordinate charts. (This is the number of permutaions that actually 
change e). Again, since the co-ordinate transformations are analytic functions, this shows 
that Grm,M is a complex manifold of dimension m{M — m). 

2. 3. The Symplectic Structure on the Grassmannian 

We will be interested in physical systems for which the Grassmannian is a Phase space. 
Since it is compact, it cannot be the cotangent bundle of any cofiguration space. It is not 



possibble to decompose the dynamical variables into configuration space and momentum 
space variables. A canonical formalism is still possible if we have a symplectic form;i.e., 
a closed, nondegenerate 2-form. Since each connected component is a homogenous space, 
it is natural to look for one that is invariant under the action of U{M). In fact there is a 
unique (upto overall constant) homogenous symplectic form. In terms our decsription of 
the Grassmannian in terms of hermitean matrices, 

u = --tr$(i$<i$. (25) 

The normalization is chosen so that the integral of iv over a two-sphere embedded in the 
Grassmannian is an integer multiple of 2n. This will be convenient later. 

The form u is obviously invariant under the transformation $ — * g^g' where g is a 
constant unitary matrix. To see that uj is closed, note that 

du = --tr(rf$)3 = --tr(rf$)3$2 ^26) 

8 8 

since $^ = 1. By differentiating this constraint, 

$rf$ + rf$$ = 0. (27) 



Now, 



div = --tr$(rf$)^$ = -tr(rf$)^$^ = -div (28) 



so that it is zero. Since u is homogenous, it is enough to verify nondegeneracy at one point 
in each connected component, say $ = e G Grm,M- A tangent vector U at this point is a 
hermitean matrix satisfying 

[e,U]+ = 0; (29) 



i.e., of the form 



u-(j o)- (^») 



Then, 



« « ^ t t 

I.I it — 1) 1 1 



uj{U,V) = --tYe[U,V] = --tv[u^v-v^u]. (31) 
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If uj{U, V) = for all V implies that w = 0; a; is nondegenerate. Thus a; is a sympletic 
form on Grm,M- 

In fact what we have is a special case of the coadjoint orbit construction of homogenous 
symplectic manifolds due to Kirillov. The space of hemrmitean matrices is the Lie algebra 
of the unitary group. Since the trace is an invariant inner product, we can identify vector 
space with its dual. The adjoint (identified with coadjoint) action is $ ^ 9^9 ■ The 
Grassmannian Grm,M is just the orbit of e. 

This suggests a generalization of our construction to more general coset spaces, 'flag 
manifolds'. Let ^ be some hermitean matrix with eigenvalues x^ with degeneracy mo, for 
a = 1, ■ ■ -d. Clearly J2a=i ^« = M. The orbit of ^ is the set of all elements related to it 
by an action of U{M): 

0^ = {<^ = 9C9'^\GeU{M)}. (32) 

This is a connected component in the flag manifold 

d 
F/(M) = {$|$t = $; JJ($_a;^) =0}. (33) 

q:=1 

The particular connected component containing ^ is picked out by putting enough con- 
straints on the traces to determine the multiplicities: 

d 
0^ = {$|$t = $. JJ($_^^).tr$ = ^m«a:«;---;tr$^-i = 5^m«4-i} (34) 

Qf=l 

It is enough to have d — 1 trace constraints since X]"^a = ^- ^s a homogenous space, 
O^ = U{M)/U{mi) X ■ ■ ■ X U{md)- The homogenous symplectic form defined by Kirillov's 
procedure is just 

u; = --tr$(rf$)2. (35) 
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It is clearly invariant under the action of U{M). At the point ^, 

u{U,V) = -UYi[U,V] (36) 

which agrees with Kirillov's definition [14] upto a constant. It should be interesting to 
consider the generalization of the hadron theory whose phase space is an infinite dimen- 
sional flag manifold; however, it would no longer be equivalent to 2DQCD for fermionic 
matter fields. Perhaps there is some more general quantum field theory that is related to 
it. 

2. 4- Poisson Brackets 

We recall some basic facts of classical mechanics to define the terminology. [15]. 
A symplectic form uj has an inverse uj~^ which is an antisymmetric contravariant tensor. 
Acting on a pair of one- - forms, it will produce a function. In terms of coordinates, 

UJ = (jj^j^dx^dx^ and (jJ~^ = '^^'^ gfir gfr , ^^'^ being the inverse of the matrix uj^^^. 

Observables of a classical dynamical system are smooth functions on a symplectic 
manifold, its phase space. Given a pair of functions /i,/2, their Poisson Bracket is the 
function defined by 

{hj2} = -u^-\dh,df2). (37) 

This bracket is antisymmetric and satisfies the Jacobi identity if u is closed. The Poisson 
algebra of a complete set of observables defines the symplectic structure uniquely. 

To each smooth function / their corresponds a vector field, defined by 

uj{Vf, .) = df. (38) 

Vf can be viewed as the infinitesimal canonical transformation generated by /. Then, 

{/i,/2}= CvJi=oj{Vf,,Vf,). (39) 

11 



A classical dynamical system is defined by a manifold , a symplectic form uj and a function 
H, the hamiltonian. The time evolution of the system is given by the integral curves of 
the vector field Vh- Functions that have zero Poisson bracckets with H are constant under 
time evolution. 

In our case the matrix elements of $ form a complete set of observables. Of course 
they are not all independent, being related by the constraints on $. It is convenient to 
think of the symplectic structure in terms of the Poisson algebra of these functions. For a 
constant hermitean matrix u define the function 

fu = -U^^u. (40) 

We can show that 

{fu-,fv} = f-i[u,v]- (41) 

Let us first find the symplectic vector field Vu associated to fu- Let us think of the vector 
field as a matrix valued function 14 (^) satisfying $F^^($) + 14($)$ = 0. Then, 

u;(K, .) = dU => tr (-^^ $[K($), d<^] = -^trd^u (42) 

That is, 

[$,K($)] = -4m. (43) 

Using the identity [$, [$, K($)]] = 4K($), we get 

K--zK$]. (44) 

This is just the infinitesimal action of U{M) on Gr{M). Then, 

{fujv}m = - cvjvm = /-.[.,„]($). (45) 

Now, let e*- be the Weyl matrices, (e*)^^ = 5^5^-. Then, the matrix elements of $ are given 
by $* = tr$e*-. Thus the Poisson algebra of the matrix elements can be written as 

{$},$f} = -2z($j5;-$i5^^). (46) 
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The Poisson algebra of a dynamical system whose phase space is the Grassmannian can 
be defined by the above relations on its generators along with the constraints ^l^l. = 5|.. 

The functions fu generate the infinitesimal action of the Unitary group; they are 
the moment maps on the coadjoint orbit. The dynamical system with one of these as 
hamiltonian (energy) is not very interesting, since it corresponds upon quantization to a 
free fermion system (see below). The hamiltonians of interest are quadratic functions of 
the $: 

if ($) = [-tr/i$ + ^tr$G'($)] (47) 

(The constant factor is chosen for later convenience). Here h is some constant hermitean 
matrix. G is a positive symmetric linear operator on the vector space of Hermitean ma- 
trices. In terms of components, 

((?$);. = Gfi^i (48) 

with 

i TDh 

connected component (depends only on the trace of $). Similarly, if G*^ = Cldj + 5\C^ 
it will also be a constant. Therefore the tensor G can be chosen to be traceless in all 
the indices at the cost of changing the hamiltonian by a constant, which will not affect 
the equations of motion, {h can also be chosen traceless). The pair of irreducible tensors 
h and G defines a 'quadratic function' on the Grassmannian. Quadratic functions on 
Grassmannians have been studied before, [16] but the ones we are interested in do not 
seem to have a natural geometrical meaning. Any G can be expanded in terms of its 
eigenvectors, G^^^ = Yla'^'^^ak^Vi^ 7" being positive numbers. The particular tensors G 
that arise in our theory have the Weyl matrices as eigenvectors. It would be interesting to 
obtain an algebraic or charecterization of these tensors, independent of the connection to 
fermionic systems. 
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Gfi = Gfj = G*^\k (49) 

and Glh^*-' ^i > 0. If G'w = B)5f^ +6)BK the quadratic function will be a constant on each 



The extrema of the Hamiltonian function (static solutions of the equations of motion) 
on the Grassmannians are of special interest. A variation that preserves the constraint 
$^ = 1 is of the form 5$ = [$, U], where U is an arbitrary anti- hermitean matrix. Thus 
the condition for an extremum is 

[/i -(?($),$] =0. (50) 

If (j = 0, the solutions are of the form $ = Xli ^i^i ® ^ ' with e^ = ±1 and Ui being 
the eigenvectors of h. In each connected component, the minimum is corresponds to the 
choice Ci = — 1(+1) for the lowest m eigenvalues of h. Even the case G = corresponds 
to a highly nonlinear classical system, since the phase space is a curved manifold. The 
minimum in Grm,M has the physical meanining of the ground state of a free fermionic 
system, with m fermions. The operator A^ = |(1 — $) (which satisfies A^^ = A^) has as its 
eigenvalues the occupation numbers of the fermionic system. The additional (quadratic) 
term represents interactions between the fermions; this can of course, distort the ground 
state. The configuration $^ with the minimum energy in each connected component 
Grm,M is the semi-classical approximation to the ground state of the fermionic system in 
the expansion we will introduce soon. 

The time evolution of our system will be determined by the ordinary differential 
equations 

If = [ft -6(4).*]. (Bl) 

Small oscillations around the minimum $o in Grm,M are of physical interest; they corre- 
spond to charge density waves (mesons) of the fermionic system. (A small fluctuation will 
grow exponetially in time if we expand around a critical point that is not a minimum.) 
The equation for small fluctuations is 

l^ = [h- G{^o).m + [$o, G{6^)]. (52) 

2 at 
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In other words, the smaU osciUations satisfy ^^ = i^((5$), K being the hnear operator on 
the tangent space of the Grassmannian at $0? 

k= ad (/i - ^($0)) + ad ^qG. (53) 

2. 5. Action of the Grassmmanian system 

By a straightforward apphcation of the usual point of view of on multivalued actions, 
we can construct an action principle that defines the classical system on the Grassmannian 
completely. 

In hamiltonian mechanics, the symplectic form is usually expressed in canonical coor- 
dinates, u) = dq^ Adpi. Given a path in phase space c(t) = {q^{t),pi{t), its action is defined 
to be 

S= fp.dq'- f H{p{t),q{t))dt (54) 



the subscript indicating that integrals are evaluated along the path c. The path along 
which this is stationary is the classical trajectory. If the curve c is closed, this can also be 
expressed as 

S= I u- I H{p{t),q{t))dt (55) 

J s J c=ds 

where s is a surface in the phase space whose boundary is c. (We will be interested only 
in situations where such a curve exists) . The latter point of view is suited to situations 
where the symplectic form is closed but not exact. 

This gives us an action principle on the Grassmannian defining our system: 

S = -- tr^{s,t)----dtAds- / tr[-/i$(t) + -$(t)G($(t))]rft. (56) 

Here $(t) is a periodic function of t decsribing a closed curve on the Grassmannian. Each 
component of the Grassmannian is simply connected, and $(s,t) is a deformation of this 
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curve to a point (i.e., it is a surface whose boundary is $(t)). Since ui is closed, a continuous 
deformation of the surface itself will not change the value of the action. However, if it is 
changed by a closed surface, the value of the action can jump by a constant. Since u has 
been normalized such that its integral over any closed 2-surface is an integer multiple of 
27r, the change of S will also be an integer multiple of 27r. This has no effect on the classical 
theory, but does affect the quantum theory. In order that the quantum path integral 

V^e-^^ (57) 

have a single valued integrand, -^ must be an integer. This is the integer denoted by N^. 
below. We will see that the same restriction arises in the canonical quantization method. 

3. Quantization on the Grassmannian 

We have been studying a classical system whose phase space is the Grassmannian 
with the hamiltonian 

if = tr[-/i$ + ^$G'($)]. (58) 

We will now show that the quantization of this system yields a system of fermions with 
a hamiltonian that is quartic in the fermion operators. In this sense, the Grassmannian 
system we described earlier is the classical limit of a fermionic system. The classical limit 
of the fermionic system corresponds to the one in which the occupation numbers go to 
infinity. 

There are several different routes to constructing the quantum theory. The simplest 
is the algebraic (or canonical) quantization. 

In the algebraic point of view, the classical system is defined by the (i) Poisson brackets 
of a complete set of observables 

{$},$f} = ($i5^-$)5n, (59) 
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(ii) the constraints satisfied by them, 

^H = ^l (60) 

and (iii) the hamiltonian 

i7 = tr[-/i$ + -$(?($)]. (61) 

One way to pass to the quantum theory is to associate to each classical observable / 
(function on the phase space) a self- adjoint operator / on a quantum Hilbert space such 
that 

{h.h\ = in{hj2\. (62) 

In conventional canonical quantization, one thinks of the canonical coordinates as the basic 
observables and require them to satisfy the Heisenberg relations. For more complicated 
operators, it is possible to preserve the connection of Poisson brackets to commutators 
only upto higher order terms in Ti. The real parameter li measures how much the quantum 
theory deviates from the classical theory. 

It is not a good idea to express our theory in terms of canonical coordinates, as 
such coordinates would be singular somewhere in the phase space. Instead, we will take 
the Poisson brackets of the $* as the analogues of the canonical commutation relations. 
Quantization then consists of finding a representation of the above algebra in terms of 
operators $* on a complex Hilbert space satisfying 

\^)M\ = -inm)-^b\). (63) 

This means that the operators -^ provide a representation of the unitary Lie algebra. We 
can assume this representation is irreducible since otherwise, the whole theory will split up 
into a sum of irreducible components that do not interact with each other. The real and 
imaginary parts of $*• are observables, so that the corresponding operators must satisfy 

17 



the hermiticity condition 



$]' = M (64) 



analogous to the classical condition $** = $*. This condition requires the representation 
of the Lie algebra to be unitary. 

Unitary irreducible representations of the Lie algebra U{M) are well known. A basis 

in any irreducible representation is given by the simultaneous eigenvectors of ^$^ ( no sum 

on i), 

-^l\w >= Wi\w > . (65) 

n 

In a unitary representation, these weights Wi are integers. There is a unique lowest weight 

vector, satisfying 

^\w>=0 for i > j (66) 

Also the lowest weight is a nondecreasing sequence wi < W2- ■ -wm- There is in fact a 
one-one correspondence between such sequences and unitary irreducible representations 

of t/(M). 

At this point our quantization scheme is ambiguous since there are many irreducible 
representations. (The usual canonical quantization is more unique since the Heisenberg 
algebra has only one unitary irreducible representation). Now recall the constraints $!^fc = 
S1;^l = M — 2m of the classical theory. If we require that these be true at least in the 
sense of an expectation value on the lowest weight state. 



< w\^'j\w >< w\^l\w >= si; < w\^i\w >= M - 2m (67) 

we will get 

h^wf = l n^w^ = M-2m. (68) 

i 

This implies first of all that ?i is the inverse of an integer, ?i = -^.( This integer is 
the number of 'colors'; the reason for this terminology will become clear in a minute). 
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Furthermore, it determines the lowest weight vector in terms of this integer: 

Wi = -Nc---Wm = -Nc,Wm+l= Nc--- (69) 

Thus each connected component of the phase space corresponds to an irreducible repre- 
sentation of the unitary Lie algebra. The representation is unique upto the choice of the 
integer Nc, whose inverse has the meaning of h. In particular, the large Nc limit is the 
classical limit of the theory. 

This representation of U{M) can be decsribed more explicitly in terms of fermionic 
operators. Define a set of operators satisfying the anti-commutation relations 

K, a],]+ = S]6l K, ajt]+ = 0, [at-, a],]+ = (70) 

where i,j = 1, ■ ■ -M, and a,b = 1, ■ ■ ■ Nc- Then, there is a representation of this algebra 
on a fermionic Fock space JF of dimension 2^^". The bilinears 

^5 = ]^Ek"'4] (71) 

provide a representation of our algebra of $'s on this Fock space. However, as it stands 
this representation is reducible. For example, the operators 

Qt = -^K, 4] + ^^>''^ 4b\ (72) 

which generate a SU{Nc) algebra commute with the $*•. To get an irreducible representa- 
tion, we must look at the subspace of vectors annihilated by these operators. 

This SU{Nc) symmetry is called 'color' symmetry, since it has to do with the number 
of copies of otherwise identical fermions. It has no meaning within our theory, and it 
is natural that all our states be invariant under color. We emphasize that Q'^ does not 
generate an ordinary symmetry that commutes with just the hamiltonian; it is a gauge 
symmetry, that commutes with all the observables of our theory. 
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The fermion number operator al^a^"" also commutes with the $^ and so must be fixed 
to get an irreducible representation. If we restrict ourselves to the states of zero 'color' 
and of fermion number mNc, 

^mo = (1^ >e J'lQtli; >= 0; ^^aC^m^ >= ^1 >} (73) 

we get an irreducible representation. The lowest weight state is the state in which the 
states labelled by z = 1, • ■ -to are occupied, while the remaining ones are unoccupied 

alJw >= for z < m; aia\w >= for z > m. (74) 

The lowest weight can be calculated: 

$^|m; >= — y^[a'*",aia]|«' >= tIw' > for i < m or i > m (75) 

1\ r 



" c 



agreeing with the earlier result. Thus the color singlet fermionic states do indeed provide 
the particular irreducible representation we need. 

We can now show that the constraints are true as operator equations, in the Hilbert 
space J^rnO- 

^^ = Sl, tY^i = M-2m. (76) 

The linear constraint is trivial to see. As for the quadratic constraint, if we expand 

c 

it will be clear that every term can be written as a bilinear or as a bilinear multiplied by 
the latter terms also to bilinears. Thus we see that 



a-a^^ on the left or the right. Now we can use the condition of color invariance to reduce 



^^i = adl + b^l (77) 

20 



where a, b are real numbers that may depend on m. We can determine them them by taking 
the matrix elements of the equation in the lowest weight state. If we put i — k < m, the 
expectation value of the R.H.S. is a — 6. For the L.H.S., 



^ j<m 



1 llJ.^&U„^l|2_^< 



so that a — 6 = 1. If we consider i = k > m, we will get a + b = 1. Thus a = 1 and 6 = 0, 
proving our identity. 

Thus we have a representation of the commutation relations that satisfies the con- 
straint. To complete the construction of the quantum system, we find the hamiltonian 
operator. 



1 

—I 
2 



H = -h]^i + -G^^,^^^. (78) 



This describes a system of fermions interacting through a two body 'potential'. G^j^i is the 
scattering amplitude in Born approximation for a pair of fermions. When the interaction 
is strong, it becomes very difficult to solve the system in this language. It is simpler to 
solve it in the -^ expansion, which corresponds to the semi- classical approximation of 
the Grassmannian model. 

The quantum fiuctuations in the observables (real and imaginary parts of $*) go to 
zero in the large Nc limit. This is clear from the fact that the commutators of the $*■ 
vanish as -^ and the Heisenberg uncertainty relation 

(AA)2(AS)2>i|<^|[A,S]|^>|2 (79) 

for a pair of hermitean operators. For finite Nc, one can consider the states in J-'mO that 
minimize the uncertainty of measuring the various components of $*•. These would be the 
analogue of minimum uncertainty wave packets of quantum mechanics ('coherent states'). 
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One invariant measure [3] of this uncertainty is 

AC2 =< VI^5*'IV' > - < i^mH >< V'l^'IV' > • (80) 

The states that minimize this are the lowest weight states, i.e., staisfy $^|'i/' >= for some 
choice of basis in the Lie algebra. In fact the set of minimum uncertainty states is just the 
orbit of a particular lowest weight state under the action of U{M). The isotropy group 
of alowest weight state is U{m) x U{M — m) so that this orbit is a Grassmannian. We 
have here an embedding of the Grassmannian into the Hilbert space J-'mO- This extends 
to a holomorphic embedding of Gr{M) to the projective space V{J-'mo)- The points on 
this submanifold of J-'m,o form an overcomplete basis; i.e., any state can be written as a 
linear combination of the points on the orbit of the lowest weight state. For Nc = 1 { when 
the condition of color invariance is trivial) this is the standard Plucker embeddingof the 
Grassmannian into the projective space V{C^"^)). 

To summarize, the quantum Hilbert space of the dynamical system on the Grass- 
mannian is the color singlet sector of a fermionic system. -^ plays the role of h in this 
quantization. A hamiltonian on Gr{M) that is quadratic in $ leads to a fermionic sys- 
tem with four point (two-body) interactions. The points of the Grassmannian correspond 
to coherent states in the Hilbert space; in the large Nc limit, the quantum fluctuations 
dissappaear and we recover the classical model on the Grassmannian. 

4. Hadrons in two dimensions 

4. 1. The Grassmannian Gri 

We will first give a reasonably precise definition of the phase space of out theory, which 
will be an infinite rank Grassmannian [5]. There is a homotopy invariant (called 'virtual 
rank') that plays an important role (baryon number) in our theory. Without a precise 
definition it would not be possible to see that there is such an invariant. We dont need to 
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maintain this level of mathematical rigor for other parts of the paper, since the answers 
do ot depend critically on topological subtleties. It might be useful to consult a standard 
reference on operator ideals (e.g. [17]) before reading this section; the summary given 
below is very brief and not very pedagogical. 

Let ?i be a Hilbert space and 7i -t a pair of infinite dimesnional orthogonal subspaces 
such that Ti = Ti -® H +. Define the self-adjoint operator e which has eigenvalues ±1 
on 7i -I-. The 'restricted Unitary group' Ui is defined by 

Ui = {g\g^g = l;[e,g]el2} (81) 

X2 being the space of Hilbert-Schmidt operators. (An operator A is Hilbert-Schmidt if 
A' A is trace class: trAi A < 00). Thus the elements of Ui are unitary and do not mix 
the subspaces H ± 'too much'. This group is of interest in quantum field theory [5] 
because the bilinears of a 1 + 1 dimensional fermionic field theory form a representation 
of its Lie algebra. There, 7i is the one particle Hilbert space; Ti ± are the positive 
(negative) eigenspaces of the Dirac hamiltonian. Thus e is the sign of the first quantized 
energy operator. In general, we would split the one particle Hilbert space into H _, the 
states with energy less than the Fermi energy and 7i +, the states with more than Fermi 
energy. In nonrelativistic quantum mechanics the dimension of 7i _ is finite; in relativistic 
quantum mechanics, this space is infinite dimesnional. 

We will be interested in the infinite rank Grassmannian 

Gri = Ui{n)/U{n-)xU{n+). (82) 

U{ H ±) is the subgroup of operators that leave the negative (positive) energy states 
invariant. (Equivalently, U{ H -) x U{ H +) is the subgroup that commutes with e).It 
can be shown that this is an infinite dimensional manifold modelled on the Hilbert space 
X2( 7i _, 7i _|_) of Hilbert-Schmidt operators from 7i _ to 7i _|_. We will now define two 
dimensional hadronic theory as a classical dynamical system with Gri as the phase space. 
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As before the Grassmannian can be parametrized by operators $: 

Gn = {$|$t = $;$2 _l;[e,$] G J2}. (83) 

It is clear that any such operator can be reduced to e by a unitary transformation: $ = 
geg' . (The convergence condition on $ will ensure that g itself satisfies the convergence 
condition of Ui.) But g is ambiguous upto a right multiplication hy h E U{ H -) x 
t/( 7i _|_). Hence there is a 1-1 correspondence between $ and points of the coset space 
Gri . The point $ = e will be the 'vacuum' configuration of the theory so in fact it is more 
convenient to introduce the 'normal ordered' variable M = $ — e. Thus, we parametrize 
the Grassmannian by 

Gn = {M|Mt =M;M^ + eM + Me = 0; tr| [e, M] ^ < cx)}. (84) 

If we decompose M into a 2 x 2 matrix with respect to the splitting H = H -Q) Ti. + 
M = I I I , the convergence condition just says that /3 is Hilbert- Schmidt. The 
quadratic constraint now implies that a and 6 are trace- class. To see this write M = 
geg' — e where, g = { , I . Here 6, c are Hilbert-Schmidt and a and d are only bounded. 



c d 
Then, a — —aa' + bb' + 1 etc.Using the unitarity of g, {aa' + cc' = 1), we see that 

a = bb' + CC' which is trace class. This technical remark will be useful in defing the 

topological invariant called 'virtual rank' below. 

The group Ui has an infinite number of connnected components labelled by an integer. 
To see this, consider g = [ t I decomposed into blocks as before; a : Ti - ^ H -,b : 

H + -^ 7i_,c: H - ^ H +,d : H + ^' H +. Since g is invertible and 6, c compact, it 
follows that, a and d are Fredholm operators of opposite index. 

We digress a little to give some basic definitions and results of operator theory. There 
are many stanadard references where results are stated more precisely e.g., [17]. An 
operator is compact if there is a sequence of finite rank operators that converges to it in 
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the operator norm. Hilbert-Schmidt and trace class operators are in particular compact. 
Not all bounded operators are compact: for example the identity is not. An operator is 
Fredholm if it is invertible modulo the addition of a compact operator. The kernel of an 
operator a is the subspace of all |V' > such that alifj >= 0. The co-kernel is the set of 
\ip > such that a'\ijj >= 0. The dimension of the kernel is not a continuous function of a. 
However, the index of a which is the dimension of its kernel minus that of its cokernel is 
invariant under continuous deformation of a. 

Returning to our context, the index of the submatrix a of g is the net number of states 
in 7i _ which are pushed out to 7i _|_ by g. The index of a is invariant under a continuous 
deformation of g. The index of d is just minus that of a: as we change g, if a certain 
number of states cross from 7i _ to 7i -|-, an equal number will cross in the opposite 
direction. Two elements of Ui with the same index are connected by a continuous path: 
Ui is the disjoint union of connected components labelled by the index. 

It is easy to give an example of a unitary operator g with non- zero index for a. 
Suppose we label an orthonormal basis of 7i by integers, with Cn for n < (n > 0) 
spanning 7i _( 7i +). An element of the connected component of Ui with index k is the 
shift operator a^ 

(^ke-n = e^+fc- (85) 

Note that if we had not imposed the condition that [e, g] is Hilbert-Schmidt, the Unitary 
group would have been contractible; there would be no such topological invariant. We 
will see that this topological invariant has the meaning of baryon number in our hadronic 
theory, which is an essential feature. More generally 'n'2m{Ui) = Z for even homotopy 
groups and 7r2m+i{Ui) = for odd homotopy groups. 

The Grassmannian is also a disjoint union of connected components labelled by an 
integer. This was already true for finite rank Grassmannians. In fact — |tr($ — 1) is the 
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dimension of the subspace in which $ has eigenvalue —1 ( the rank of the Grassmannian) . 
For our present case this would diverge. We can get a convergent quantity by subtracting 
the rank of the vacuum $ = e; i.e., consider — |tr($ — e) = — |trM. This is the 'virtual 
rank' of the Grassmannian, the difference between the number of negative eigenvalues 
of $ and those of e. (A technical remark: we can define trM to be tra + tr5. Here 
M = I I I in a basis in which e = I p. ^ ) • We showed earlier that a and 5 are 
trace class. This indirect definition of trM is necessary since M itself is not trace class.) 
In fact if M = geg' — e, it is easy to see that — ^trM = index (a). (For example, if g = ai 
the shift operator, (e^, aieai — ee^) = sgn (n — 1) — sgn (n). Then trM = —2 and the 
index of a is 1.) 

It is best to think of TC as the space L'^{R;C) of complex valued functions on the 
real line. In this parametrization, e is the sign of the momentum operator. 



Here, 



hdky'''' sgn (k). (86) 



[dk] = ^. (87) 

Ztv 



We can also describe it by its integral kernel in position space 

i / 1 

e{x,y) = —V 



27r \x — y J ^ 

This is, i times the well-known Hilbert transform operator. The Grassmannian can be 
parametrized by the integral kernel M{x,y) of the operator M. The condition of self- 
adjointness is just 

M*ix,y) = M{y,x). (89) 

The quadratic constraint is the nonlinear integral equation 

dy[M{x, y)M{y, z) + e(x, y)M{y, z) + M(x, y)t{x), z)\ = 0. (90) 
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The convergence condition is best understood in terms of the Fourier transform 

M{p,q)= f e'P''-''^yM{x,y)dxdy. (91) 

The condition is that the off-diagonal component of M with p and q having opposite sign 
has finite norm: 

/ [dp] [ [dq]\M{p,q)\^ <oo. (92) 

Jp>0 Jq<0 

The topological invariant ('virtual rank') of the Grassmannian is then 

w{M) = -^ / dxM{x,x) (93) 

which we will be shown to have the physical meaning baryon number. 

4. 2. Symplectic structure and Poisson Brackets 

We saw in the finite rank case, that the Grassmannian is a symplectic manifold with 
u = — |tr$(rf$)^. The analogue in infinite dimensions is just uj = — |tr(M + e){dM)'^. A 
tangent vector at the point M = is a matrix of the form U = i + J with tru^u < 00. 
It is straightforward to see that uj{U, V) = —^tr{u'v — v'u). This shows that u exists and 
is non-degenerate at the vacuum. That u exists and is non-degenerate at an arbitrary 
point follows from its homogenity (see below). The proof that du = goes through as 
before. 

The transformation $ -^ 9^9' or M ^ 9Mg' + 9[e:9'] leaves the symplectic form 
invariant. From general principles, there must be functions that generate the corresponding 
infinitesimal canonical transformations, y^t = i[uje + M]. (The matrix u = i ,j , J must 
have diagonal elementsa, b bounded and off-diagonal elements c, (i G X2 in order that e*^ 
be in Ui.) That is, there must be a function f^ such that uj{Vu, ■) = dfu- In the finite 
dimensional case, this is just fu{^) = trt(,$. In the finite dimensional case, these are just 
linear functions of $. For our present case, again tr-u^ will diverge; we must subtract 
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the vacuum contribuition and consider instead fu{M) = —^truM = — itrw($ — e). (It is 
easy to check that the trace in fu{M) does indeed converge.) But we wiU pay a price for 
this vacuum subtraction. The Poisson brackets of the fu will no longer be that of the Lie 
algebra of Ui: there will be a non-trivial central extension. 

The computation that verifies that fu{M) generates Vu is straightforward: 
^(K, .) = -^tr(e + M) [Vu{M), dM] = ^tr[K, e + M]dM 

o o 

= -tr[e + M, [e + M, u]]dM = --trudM = dfu- 
8 2 

Now we get 

{fuJv}{M) = - CvJv{M) = -itr[u,e + M]v = /(_,)[,,„](M) + ^tr[e,u]v. (94) 

This is the central extension of the Lie algebra of Ui . We cannot remove the central term 
by redefining fu{M) -^ fu{M) + treu since the second term does not converge in general. 
This nontrivial extension of the unitary Lie algebra was found in a different context by 
Kac and Petersen. In the co-ordinate basis, these Poisson brackets are, 

^{M(x, y), M{z, u)} = 6{y - z)[M{x, u) + e{x - u)] - d{x - u)[M{z, y) + e{z - y))]. (95) 

This may be verified by multiplying both sides by u{y,x)v{u,z) and integrating over 
X, y, z, u to get the ealier Poisson bracket relation. 

Occasionally the above Lie algebra is also called the Woo algebra. This algebra arises 
in the matrix model approach to string theory. Our theory of hadrons is a nonlinear 
matrix model, the dynamical variable being M{x, y). We believe that in fact our theory is 
a string field theory in two target space dimensions written in light cone gauge. Finding an 
ungauge fixed (manifestly reparametrization invariant form) of this theory is an extremely 
interesting problem. It is already known that the linear approximation to our theory 
(which is described by 't Hooft's integral equation) is equivalent to the free Nambu string 
in two target space dimensions [18] 
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The Poisson brackets above along with the constraints can be thought of as the def- 
inition of our phase space. The interpretation in terms of the Grassmannian shows how 
natural these Poisson brackets are, given the constraints. Conversely, given the above Pois- 
son brackets, we can look for homogenous symplectic manifolds on which the Lie group 
acts. These are given by the co-adjoint orbits, of which the Grassmannian is one of the 
simplest. The other co-adjoint orbits do not describe fermion fields upon quantization. 

4. 3. The Hamiltonian 

We will now describe two dimensional classical hadron theory. The most convenient 
co-ordinate system of Minkowsky space for our purpose is the one with metric ds"^ = 
du{du + 2dx). In terms of conventional space- -time co-ordinates ds^ = (dx^)^ — (dx^)^ 
and u = x^ — x^,x = x^. The surface u = constant is a null line and x is a co- ordinate 
in this null direction. The vector ^ (for fixed x) is just the same as the time-like vector 
-Af. Thus momentum component Pu is just energy while the x component of momentum 
(which we will call just p) is p = po + pi. The mass shell condition for a particle of mass 

fx is 

1 A*^ 

2puP - P^ = l^'^ ^ Pu = -[p ^ ]• (96) 

/ p 

The Lorentz transformation is, in this co-ordinate system, 

u ^ e u X -^ ~ sinh 9u + e~ x (97) 

where 6 is the rapidity. From the invariance of p^du + pdx, we get 

p ^ e p Pu ^>- e~ Pu + sinhOp. (98) 

The phase space will be given by the initial conditions on the surface u = 0; the dynamics 
will be given by the hamiltonian function H which will generate translation in u. There 
should also be functions L and P generating Lorentz tranformations and translations, 
satisfying the Poincare algebra: 

[P, L] = [L, P]=P [L, H] = -H + P. (99) 
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A dynamical system is completely specified by the phase space (a manifold along with 
a symplectic form) and a Hamiltonian. We will define two dimensional hadron theory by 
choosing {Gri,uj) as the phase space. Thus our dynamical variable is a function of two 
points M{xj y) satisfying the constraint 

M*{x, y) = M{y, x), j dy[M{x, y)M{y, ^) +e(x, y)M{y, z)+M{x, y)e{y, z)] = 0. (100) 

This specifies the initial condition at m = 0. The value at any later vaue of u is given by 
solving the classical equations of motion 

dM{x,y) 
du 

We choose the hamiltonian to be 



{H,M{x,y)]. (101) 



H = I dxdyh{x - y)M{x, y) - -f I dxdyG{x - y)M{x, y)M{y, x). (102) 

The operators h and G are defined by 



-{p + ^y, G = - 

2 p p^ 



^=nb+— ); G=-, (103) 



where p = —i-j^- The theory depends on the parameters g and (i with the dimensions of 
mass. We will show later that n is related to the quark mass and g to the gauge coupling 
constant of 2DQCD. The integral kernels of the operator are 

I , i 1 

h{x - y) = -[-i5' {x - y) + - sgn {x - y)], G{x - y) = --\x\. (104) 

The choice of h as the kinetic energy operator is motivated by the earlier discussion of the 
mass shell condition. 

Although the field variable is bilocal, the theory we have just defined is a Poincare 
invariant field theory. Translation invariance in the x direction is obvious. Under Lorentz 
tranformations, the field transforms as 



^1 



M{x,y)^e^M{e^x,e%). (105) 
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It is straightforward to check that H and P = trpM transform as the components of a 
Lorentz vector under this transformation. Once we assume the form () for the hamiltonian, 
specific expressions for h and G are required by Lorentz invariance alone. It should be 
possible to show by an argument based on string theory that () is the only possible form 
for the hamiltonian. 

The equations of motion that follow from this hamiltonian are, 



2 du 2 

+ ^M dz[G{y - z)e{x, z)M{z, y) - G{z - x)e{z, y)M{x, z)] 



+ ,^/..M(.,.)M(.,.)[G(.-.)-G(.-.)l. 
By construction, the quadratic constraints are preserved under time evolution. 

The solution M = is the vacuum solution. Small oscillations around this are de- 
scribed by the linear approximation to the above equation. These will describe mesons. 
There are also stationary solutions whaich are very far from this vacuum solution, which 
describe baryons. The quantity B = — ^trM, is an integer valued topological invariant and 
hence is conserved by time evolution. The lowest energy solution with B = 1 should be 
identified with the baryon. Linearizing (105) around this stationary solution will describe 
meson-baryon scattering; the bound states in this linearized equation will dsecribe excited 
states of the baryon ( analogues of A, A^* etc.). 

4- 4- Linear Approximation 

Let us now study the linearization of the above theory around the vacuum M = 0. 
We must linearize not only the equations of motion but also the constraints. If we drop 
the second order term, the constraint becomes, 

J dy[eix, y)M{y, z) + M(x, y)e{y, z)] = 0. (106) 
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Since our equations are translation invariant, it is better to use the Fourier transformed 
variable 

M{p,q)= f dxdye'^-P''+''y^M{x,y). (107) 

The constraint becomes: 

[ sgn (p) + sgn {q)]M{p, q) = 0. (108) 

In addition the hermiticity consition becomes 

M*{p,q)=M{q,p). (109) 

This means that our field operator has only off-diagonal components; i.e., those with 
opposite signs for p and q. (Of course, this is just the statement that M is tangential to 
Gri at the origin). The translation, M{x, y) -^ M{x + a,y + a), and 

M{p, q) -^ e^(-«"+P"^M(p, q). (110) 

Thus P = p — q has the meaning of total momentum. Due to the hermiticity condition, it 
is sufficient to consider the case p > and g < 0, so that the total momentum is positive. 
Define now the dimensionless variable ^ = p/P; then ^ varies in the range < ^ < 1. We 
can use P and ^ as our independent variables. Set 

x{P,x)=PM{CP,{l-OP). (Ill) 

The linearized equations of motion is, in position space, 

. dM{x,y) ^ I ^^^^^^ _ ^^^^^^ ^^ _ ^^^^ ^^^^^^ ^^^ 

+ f I dz[G{y - z)t{x, z)M{z, y) - G{z - x)t{z, y)M{x, z)]. 
In momentum space this becomes, 

i gl' = IHp) - h{q)]M{p, q)+f[ sgn (p) - sgn (q)] J ^M(p -r,q- r). (112) 
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It is clearly natural to assume that the u dependence is exponential 

M{p,q,u)^e-'^-''M{p,q) 
. In terms of the variable x(P, ^) we get the eigenvalue equation 

where M^ = 2PuP — P^ is the mass^ of the meson. This is precisely 't Hooft's integral 
equation. He has shown that all solutions are bound states ( discrete eigenvalues) and that 
spectrum is assymptotically M^ ~ ng^. 

4. 5. The baryon solution 

If ours is to be a complete theory of hadrons, it must contain baryons as well as 
mesons. From Witten's [8] arguments one should expect baryons to arise in a Hartree 
approximation. On the other hand, in the Skyrme model, they arise as solitons in a 
theory whose small oscillations are mesons. Normally, this meson theory is only known 
approximately. However, we have an exact theory of hadrons in two dimensions. This 
gives us an oppurtunity to study the baryon within this exact theory. Also, we will get a 
description of a topological soliton in a bilaocal field theory. Since we should expect our 
bilocal theory to be a string field theory, we will also get a glimse of how a topological 
soliton looks like in string theory. 

We already showed that there is a topologically conserved quantum number w{M) = 
— |trM in our theory. It takes integer values. (It will be shown later that this indeed 
corresponds to baryon number of 2DQCD.) We will look for a configuration that minimizes 
the energy, and with w{M) = 1. Unlike the vacuum solution, it will have non-zero energy: 
it is a topological soliton. Upon quantization, its mass will seen to be 0{Nc), since -^ has 
the meaning of h in the quantum theory. We will see that the equation satisfied by soliton 
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field configuration has a natural interpretation in terms of a Hartree approximation for the 
quark wavefunction. 

Our first task is to produce an ansatz that has virtual rank one. It will be much more 
convenient to deal with functions of one variable rather two, so we will seek a separable 
ansatz of the form 

M^{x,y) = -2^{x)r{y)- (114) 

In order to have virtual rank one we must require the function if) to be normalized: 



dx\il){x)\'^ = 1. (115) 

Now let us see the condition on •0 implied by the quadratic constraint, M^+eMi+Mie = 0. 
It is clear that il){x)il)*{y) is the kernel of a projection operator, so that Mf = —2Mi. Thus 
we get eMi + Mie = 2Mi. In terms of t/j this is just the condition et/j = tj). That is, tj) is 
an element of 7i _(-, the space with eigenvalue +1 for e: 

e{x-y)ij{y)dy = ijj{x). (116) 

Thus our ansatz Mi is indeed an element of Gri with virtual rank one.( It is trivial to 
check that [e, M] is a rank one operator, so that the convergence condition tr| [e, M] p < cxd 
is satisfied). 

In fact Ml is unchanged under the change of phase ijj{x) -^ e^^ip{x). The ansatz gives 
a map from the projective space V{ Ti +) to Gri. We will see that 7i _ is the set of 
negative energy states for the quarks, which are occupied in the ground state of 2DQCD. 
To produce a baryon, we must put Nc quarks in one positive energy state as well. In the 
ground state, all these quarks will have the same wavefunction, apart from a color factor 
which is completely anti-symmetric. This is the meaning of if) in the quark model. Note 
that the condition that •0 is a positive energy state arises naturally in the above formalism. 
It will be crucial for the stability of the soliton. 
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The energy of this configuration is obtained by substituiting the ansatz into the hamil- 
tonian (). The first term in the hamiltonian becomes trhM = {t(j,ht(j) or, J'i(j*{x)h{x — 
y)'ilj{y)dxdy. Noting also that M{x^y)M{y,x) = \'ilj{x)\'^\'i/j{y)\'^ we get 

if (V) = jr {x)h{x - y)^{y)dxdy - ]^f j G{x - y)\ij{x)\^{y)\^dxdy. (117) 

It is now clear that the first term is just the kinetic energy of the quark wavefunction. The 
second term is the potential energy due to a Coulomb field generated by the quarks. There 
is an attractive linear potential between the quarks. This is exactly what we would have 
obtained in a Hartree approximation. There are N^. quarks in a baryon, whose color indices 
must anti-symmetrized to get a singlet state. The quarks being fermions, the wavefunction 
must be symmetric in the remaining ( position and fiavor) indices. Thus effectively they 
are bosons and the ground state wavefunction of all the quarks will be the same. The 
second term in the energy then describes the interaction with the mean field. 

, 2 

The condition eif) = if) ensures that the kinetic energy (-0, hif)) is positive, h = ^ [p+ —] 
is not in general positive; 7i + is the susbspace in which it is positive. If it had not been 
for the constraint that ifj E 7Y +, the energy would not be bounded below and our soliton 
would have been unstable. The constraint is best understood in momentum space. If we 
define the Fourier transform 

«p)^/^(.)e-.. (118) 

we see that it has support only for positive p. In terms of position space, 'iIj{x) must be a 
function that can be analytically continued to the upper half of the complex plane. 

Thus we must minimize the energy over all configurations i^{p) satisfying 

ip{p) = for p < / \ijj{p)\'^[dp] = 1. (119) 

Jo 

The Kinetic energy is simple to understand in momentum space while the potential energy 
is simpler in position space. The energy is 

H{ip) = - / [p + —]\'>P{p)\^[dp\ + -f / V{x)\iP{x)\^dx. (120) 



2 Jo " p''^'"' '^' 2- 
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Here V{x) is the solution to the diffrential equation 

V"{x) = \ij{x)f (121) 

with the boundary conditions V^(0) = and V'{0) = 0. This variational problem can be 
reduced to a nonlinear integral eigenvalue problem. We have not been able to solve this 
problem analytically. The problem was solved by numerically minimizing the energy in 
Ref. [11]. 

Here we will find the qualitative behaviour of the soliton using a simple variational 
approximation. An ansatz that satisfies the positivity of energy and the normalization 
condition is V'(p) = Npe~^"', where a > is a variational parameter. (The even simpler 
possibility V'(p) = Ne~P"' has to discarded because it has infinite potential energy). Then, 
ip{x) = -^ (x+ia)'^ • Note that this has a double pole on the lower half plane, but is analytic 



27r (x+ia)2 

271- 



in the upper half plane. The normalization condition is -g— = 4a^. The kinetic energy 



integral is easily done: 



ll\p+^mpndp] = l[^+f^'a]. (122) 



The differential equation 



has solution 



2ri^ 1 

V'ix) = mx)f = — (123) 



1 X 

F(x) = — xarctan — . (124) 

TV a 



Note that V{x) ~ ^\x\ as \x\ -^ oo, as expected. Thus the potential energy becomes 

^f fv{xM{x)fdx = hfa (125) 



where 

4 [°^ xarctanx 



TT Jo 1 + X^ 
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dx = 2. (126) 



Thus we estimate the energy of the sohton to be 

Hi'^-\\h'^\ (127) 

where the quark mass /i is replaced by the effective mass: 

lils = li^ + 2f. (128) 

As the mass of the quark goes to zero, the hghtest meson also has a very small mass [7] . 
But we see that even as the mass of the quark goes to zero, the mass of the baryon remains 
finite! The effective mass of the quark in the baryon is heavier by an amount proportional 
to the gauge coupling constant. That is, the 'constituent quark mass' is of the order of 
the gauge coupling constant, even as the current quark mass goes to zero. We have just 
established a phenomenon analogous to chiral symmetry breaking in two dimensions. (Of 
course the chiral symmetry itself cannot break spontaneously in two dimensions; the meson 
of small mass also has small coupling, so that it is not a Goldstone boson.) 

In the quantum theory, there is an overall factor of N^. = -^ in front of the action ( 
and hence the hamiltonian) so that the mass of the baryon is N^. times the minimum of 
H as determined above. Thus our estimate /Ue// is the mass of the baryon divided by Nc, 
which can be thought of as the constituent quark mass. In the case of mesons, the energy 
of the vacuum configurations are zero so that the masses of the mesons come from small 
oscillations around the vacuum. These are therefore h = ^ order smaller than the baryon 
mass, hence are 0{N^). 

5. Quantization of Hadron Theory 

As in the finite dimensional case, we can find an action principle for QHD. The 
quantum path integral would be 

V^e~^^ (129) 
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where 



f dM dM r 

S= oj(-—, -—dt Ads- H{M{t))dt. (130) 



Again, M is an extension of M to a function of two variables t, s such that the boundary 
value is the closed curve M{t). As before, H'^{Gri) = Z and the integral of a; on a 
generator of H'^{Gri) is 27r. Thus ^ = Nc must be an integer in order that e~^^ be single 
valued. Since at the moment it is difficult to define such path integrals rigorously, we will 
follow instead a canonical (algebraic) point of view. The definition of the path integral 
can be accomplished by expoliting localization formulae. 

We can now quantize the hadron theory by algebraic (canonical) methods. We would 
convert the Poisson brackets into commutation relations: 

[fuJv] = ih{f-i[u,v] +tr[e,u]v). (131) 

Equivalently, we look for operator-valued distribuitions M(x, y) satisfying, 

[M{x,y),M{z,u)] = -ih{d{y- z)[M{x,u) + e{x,u)]- 

6{u-x)[M{z,y) + e{z,y)]). 
These will provide a a representation for the Lie algebra of Ui on some space J-'. (Each 
matrix element M{x, y) of M must itself be an operator on JF.) 

The representation we pick must be unitary: 

M{x,y)'f =M{y,x) (132) 

and irreducible. (If the M{x,y) are to be a complete set of observables, the only operators 
that commute with them must be multiples of the identity; but then the representation 
is irreducible.) Also the quadratic constraint must be satisfied at least upto terms that 
vanish as ^ ^ 0. 

Furthermore, the hamiltonian must become a well-defined (self-adjoint) operator that 
is bounded below. Even for the simplest case of the hamiltonian with ^ = this imposes 
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a nontrivial constraint on the representation. It is best to discuss this case in momen- 
tum picture, so we define M{p,q) = J dxdyM{x,y)e~'^P^~^'^y. The commutation relations 

become, 

[M{p, q) , M{r, s)] = -iTi [2'K5{q - r) [M(p, s) + sgn {j))6{p - q)] 

— 2n5{s — p) [M{r, q) + sgn {r)6{r — s)]] . 

It is easy to recognize this as a representation of the central extension of of the unitary 

Lie algebra, if we write this in terms of -^ : 

[ ^ , ^- ] = -« [2Tid[q - r) [ + - sgn {j))d[p - q)\ 

- 2tx5{s - p) [MhA + 1 sgn {r)5{r - s)]] . 
Note that the central terms are proportional to ^. 

The hamiltonian is ^o = f[dp]h{p)M (p, p) In the case of interest to us, the operator 
h has an infinite number of eigenvectors with a negative eigenvalue: elements of 7i -. It is 
clear that acting with M{p, q) on any state will add p — q to its momentum. By operating 
with M(p, q) with p < q repeatedly, we may be able to construct states of arbitarily 
negative momentum (and hence energy). The only way to avoid this catastrophe is to 
eventually arrive at a state ('vacuum') |0 >G JF such that 

M{p,q)\0>=0. for p<q (133) 

It is clear that M(p, q) forp<q can be thought of as the negative roots of the Lie 
algebra of Ui: we have just shown that the representation of JF must be a lowest weight 
representation. If the representation is to be irreducible, the lowest weight vector must 
be unique. Also, every vector in JF can be written as linear combination of the vectors 
M{pi,qi ■ ■ ■ M{pr, q'r)|0 > obtained from the vacuum by positive roots:pi > qi ■ ■ -pr > qr- 
The vacuum itself is an eigenvector of the diagonal generators: 

^li^\0 >= w{p)\p > (134) 

n 
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(To make the theory mathematicaUy rigorous, we can suppose that space is a circle so 
that the momenta take on just discrete values. At the end the radius of this circle can be 
taken to infinity.) Lowest weight unitary representations are classified by the eigenvalues 
of the diagonal elements (lowest weight) [19]. If this representation of the Lie algebra can 
be exponentiated to one of the Lie group, the weights w{p) must be integers. Moreover, 
^ is an integer A''^ for the central term to be exponentiated. (In fact the condition that 
there be a positive invariant inner product already requires w{p), Nc to be integers). The 
quantization of Nc is analogous to the quantization of the level number in the classification 
of lowest weight unitary representations of the affine Kac-Moody algebras. Indeed, we saw 
that ^ appears only as a coefficient of the central term of the Lie algebra. 

Now we are interested in representations in which the quadratic constraint on M are 
preserved in some sense. As in the finite dimensional case, we could simply impose the 
quadratic condition on the expectation value of M on the lowest weight state. This will 
lead to the condition 

[nw{p) + sgn (p)]2 = 1. (135) 

If we restrict to a topological sector where the baryon number is B we have the linear 

constraint, 

— 5Z < w\M{p,p)\w >=B ^ "^[^wip)] = B. (136) 

p p 

These conditions uniquely determine the lowest weight w{p) in terms of N^, and B: 

w{p) = Nc[ sgn {p - pf) - sgn {p)] (137) 

the 'Fermi momentum' pp being fixed to be pp = B. 

This particular representation of the unitary Lie algebra can be written in terms of 
fermionic variables, just as we did in the finite dimensional case. Introduce the fermionic 
fields satsfying the canonical anti- commutation relations: 

[r{x), xUx)]+ = s^S{x - y), [x"(x), x'{y)]+ = 0, [xl{x),xl{y)]+ = 0. (138) 
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Here a, 6 = 1 ■ ■ ■ A''^ is the 'color' index. Then, we define the lowest weight state ( vacuum 
of the free field theory) by 

r{p)\0 >= 0, for p > pF, xlip)\0 >= for p < p^, (139) 

where x is the Fourier transform of x- We can now verify that 

M(^,y) = ^:xl(^)x"(y): (140) 

satisfies the commutation relations above. The central terms arise from the fact that M 
has to be normal ordered with respect to the above lowest weight state. 

This representation of the unitary Lie algebra on the fermionic Fock space is reducible; 
in order to get an irreducible representation, we must impose the conditions 

Qa\^>=0 (141) 

on all the physical states, Q^ = /[: Xa{x)x^{x) — ■^^bXcX^i^)]dx being the 'color charge'. 
On each subspace of color singlet states of fixed baryon number, our algebra has an irre- 
ducible representation. In terms of the fermion fields, baryon number is 

1 I' ,t, 



B^Jf I :xi(i-)x°W:dl. (142) 



Now we can express the hamiltonian in terms of M and hence x- 

H = dxdyh{x - y)M{x, y) - -f I dxdyG{x - y)M{x, y)M{y, x). (143) 

As before, the operators h and G are defined by 



1. I^\ ^ 1 
-p+^; G=- 

2 p p^ 



^=n(P+— ); G=-^ (144) 



where p = — ^^- We have constructed the representation and the normal ordering of M 
such that the first term is finite. The interaction term involving G{x — y) is also now finite. 
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(If the coupling constant had been dimensionless as in the Thirring model, this would have 
been true only after a coupling constant renormalization.) 

Next we will show that the above quantization of the hadron thoery is equivalent to 
two dimensional QCD. More precisely, we will show that it is equivalent to the color singlet 
sector of QCD. As part of this we will see that A''^ corresponds to the number of colors 
and B to the baryon number of QCD. 

We will find it convenient to formulate QCD in the light- cone type co-ordinates 
defined earlier: 

ds^ = du{du -\- 2dx) , u = x^ — x^ ,x = x^. (145) 

The action of 2DQCD is 

S = ji- \F^,bFr' + qabn-<d^ - gA;,) - mS^]q'}d^x (146) 

We are considering the theory with just one flavor and SU{Nc) colors. Thus A\ are 
traceless hermitean matrices, with a,b = 1 ■ ■ ■ Nc- The generalization to several flavors is 
quite straightforward. The variables q, q are to be viewed as anti-commuting (Grassmann- 
valued) in a path integral approach. 

The Dirac matrices satisfy the relations 

(7")' = 0, (7")' = -l, [7",7l+ = 2. (147) 

A convenient representation is 

It will will be convenient label the componets of the quark flelds as follows, 
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Then, 

9 = 9V = ^(r7t x^). (150) 

As for the gauge fields, it is convenient to choose the hght-cone gauge, Aj^ = and denote 
the remaing component by A^ = A. With these choices, the action of the theory becomes, 

/-I 11 

{--tr{d^A)^+xH-idu-gA]x+-[-ivdxV+^x'd^x] + -^rn[x^r]+r]^x]}dxdu. (151) 

Now it is clear that the fields r] and A do not propagate; there equations of motion can be 
used to eliminate them in favor of x'- 

-id^V + mx = 0, dlAt + gxlx'' ^ 0. (152) 

After this elimination, the action becomes 

Here 

Ptix) = xUx)x''{x) - l-6^xkx)x%x) (154) 

is the color charge density. Also, h and G are defined as before. We already begin to see 
the sort of expressions that appear in QHD. 

The first term in the action implies that in canonical quantization, x cind x ^ire 
conjugate; the remaining terms determine the hamiltonian. The commutation relations 
are 

[rix),xUx)U = d^,dix-y), [x"(x),x'(2;)]+ = 0, [xl{x),xUy)]+ = 0. (155) 

The, Fock representation of these relations can now be constructed, based on the naive 
vacuum state 

^«(p)|0>=0, forp>0, xj(p)|0>=0 for p< 0. (156) 
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The hamiltonian is 



H = J ■ ^1^^" ■.dx+^- j G{x-y): pl{x)p''M : dxdy. (157) 

A normal ordering with respect to the above vacuum is necessary to make the hamiltonian 
of the quantum theory well-defined. It is now possible to write the hamiltonian in terms 
of the color singlet operator 

M(^,l/) = ^:xl(^)x"(?/):- (158) 

We have already seen that this operator satisfies the commutation relations of QHD. Thus 
we see the operator M(x, y) is a quark-antiquark bilinear. We see that N^. is the number 
of colors and B = J M(x, x) the baryon number of QCD. 

The kinetic energy term is already of the form tr/iM. The interaction term also can 
be written in terms of M if we use the Fierz identity [11]. This will lead to precisely the 
hamiltonian of QHD. Since the calculation has been done elsewhere [6], [11], [4], we will 
not carry it out here. Thus we see that the theory we have been studying is equivalent to 
the color singlet sector of QCD. The parameters /U, g are related to the quark mass m and 
gauge coupling constant g by the formulae, 

2 2 9 ^c ~2 2 ]\T /I ixn\ 

p =m -^^, 9 = 9 Nc. (159) 

It is to be noted that although the parameter m in the 2DQCD lagrangian is tra- 
ditionally called the quark mass, the quark is not asymptotic particle in the theory and 
has no well-defined mass. If the scattering of color singlet particles produce only color 
singlet particles, the S-matrix of QHD will be unitary. This would amount to a proof of 
confinement of QCD. It seems likely that such a proof can be made to all orders in the -^ 
expansion. 
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Appendix: Hartree-Fock Theory and Grassmannians 

1. Slater Determinants and Grassmannians 

In this appendix we will show the analogue of our approach to 2DQCD, in electronic 
physics is just the Hartree-Fock approximation. It is hoped that our point of view in terms 
of Grassmannians to Hartree-Fock theory is of some use in atomic and condensed matter 
physics. 

Let Ti = L^{R^, C^) be the one particle Hilbert space of electrons; we will think of its 
elements as complex valued functions u{x, a) of position x E R^ and spin a = ±1. * The 
Hilbert space of an m-electron system is the exterior power A"^(7i). A wavefunction in 
this space is a completely antisymmetric function t/jixi, ■ ■ ■ , Xn)- The hamiltonian of the 
TO electron system is the following operator on A"^(7i): 

m 

H = Y,[-'^l + V{xa)]+ Yl Gixa-xt). (160) 

1 l<a<b<m 

For an atom, V{x) = ^ describing the Coulomb attraction to the nucleus and G{x — y) = 
,]_ I . (We are using units with -^ = 1 and e = 1. ) The ground state is the minimum of 
the expectation value of the hamiltonian over all states of unit norm: 

Eo= inf [il^.H^). (161) 

If IIV'II = 1, the expectation value of H can be written as 

< il^\H\^ >= trhpi + trG2p2 (162) 

where h = — V^ + V and the 'one-particle density matrix' pi is a selfadjoint operator H 
with kernel 

Pl{x,y) =TO / 1p*{x,X2,- ■ ■Xm)tp{y,X2,- ■ ■ ,Xm)dX2- ■ -dXm- (163) 



* Often it will be convenient to suppress the spin variable, and let x stand for (a:, a), 
j dx for Xlcr / '^^ ^^^- The norm, for example, is then WuW^ = {u,u) = f dxu* {x)u{x) . 
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Also G2 and p2 are selfadjoint operators on the two particle Hilbert space K^iTi) with 
kernels 



P2{x,y;z,u) 



m 



t(;*{x,y,xs, ■ ■■Xni)i^{z,u,X3,- ■ ■ ,Xm)dx^- ■ -dxr 



(164) 



G2{x, y; z, u) = -[5{x - z)S{y - u) - S{x - u)5{y - z)]G{x - y). 
The one and two particle density matrices are normalized as below: 

trpi = / pi{x,x)dx = m, trp2 = / P2{x,y:,x,y)dxdy = I j. 



(165) 



(166) 



The problem of finding the ground state of this system is very hard. The idea of the 
Hartree-Fock theory is to approximate the wavefunction by a Slater determinant 

Ui{Xi) ... Ui{Xm) 



t(;{xi ■■■Xm) 



1 



^Jrn}. 



(167) 



'^ra\X\) • • • '^m\Xm) 

The one particle wavefunctions Ua ^'H are orthonormal (tta, Ub) = Sab- A change Ua{x) -^ 
9abUb{x) will change the Slater determinant only by a phase detg, so it doesnt change the 
physical state represented by it. 

The Hartree-Fock approximation to the ground state energy is the minimum over all 
choices of u of the expectation value E^^{u) of the hamiltonian in this state: 



E, 



HF 



^HF, 



- inf E"^{u). 



(168) 



Clearly, Eq < Eq^. E^^ can be woked out to be: 
E^^{u) =J2 dx dxdyul{x)[-V^ + V{x)]ua{x) + 



a=l ■ 



dxdyG{x - y) —i^[ul{x)ul{y) - ul{x)ul{y)]—-[ua{x)ub{y) - Ub{x)ua{y)]. 
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However, it is somewhat unnatural to think of the Ua{x) as the variational parameters of 
Hartree-Fock theory; a change Ua{x) -^ glub^x) by a unitary matrix does not affect the 
multiparticle state. In fact the larger m is, the larger the number of such spurious degrees 
of freedom in Ua- In quantum field theory, m is infinity, this is particularly awkward. We 
will therefore develop a new point of view in terms of Grassmannians which has a natural 
generalization to quantum field theory. Perhaps this point of view is useful even in atomic 
physics. 

To get a more intrinsic form for E^^ , let us compute pi and p2 for this ansatz. We 
will get 

pi{x,y) = P{x,y); p2{x,y; z,u) = -[P{x, z)P{y,u) - P{x,u)P{y, z)] (169) 

where 

m 

P{x,y) = J2<(^)My)- (170) 

a=l 

P is the projection operator to the m-dimensional subspace spanned by Ua- The expression 
for p2 in terms of P can be written as 

P2 = PAP (171) 

where the wedge product of two operators A, B on H is the operator on A^(7i) 

{AAB){u,v) = {Au)A{Bv). (172) 

Thus we get 

E^^{P) = trhP + trCaP A P. (173) 

In coordinate basis, 

E^'^iP) = I dx[(^ ^'2dy^^ ) + ^(^)^(^' ^)] + 

- / dxdyG{x - y) [P(x, x)P{y, y) - P(x, y)P{y, x)] 

48 



The first term represents the kinetic energy and the second term the potential energy 
due to the nucleus. The third term is the energy of Coulomb repulsion of the electrons 
('direct energy'). The last term is the 'exhange energy' due to the antisymmetry of the 
wavefunctions. (The fact that Hartree-Fock theory can be expressed in this way is well 
known [20]). The Hartree-Fock ground state energy is the minimum of E^^{P) over all 
projection operators of rank m: 

E^^ = inf E^^{P). (174) 

p2=P;trP=m 

The set of all projection of operators of finite rank defines an infinite dimensional 
analogue of the Grassmannian* 

Gr{n) = {PIP"^ ^ P-P^ ^ P]P is finite rank}. (175) 

Then trP exists and is an integer; the Grassmannian is a union of connected components 
labelled by the trace; each connected compoenent is the homogenous space 

Gr^iU) = U{n)/U{m) x U{n) (176) 

just as the finite dimensional case. Gvmi'H) is an infinite dimensional manif Each point 
in the Grassmannian represents an m dimensional subspace of 7Y, the eigenspace of P 
with eigenvalue 1. There is an embedding of GrmCH) into V^hJ^Ti), ( generalization of the 
Plucker embedding) which is precisely the meaning of the Slater determinant. To construct 
this analogue of the Plucker embedding, pick an orthonormal basis in the subspace of P: 

P^J2^a^ul (177) 

a 

and define the corresponding point in AJ^Ti to be 

Ui /\U2 /\ • • -Wjn- (178) 



* We can require P to be compact rather than finite rank; any compact projection 
operator is also finite rank. 
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Again, this vector changes by the phase det g as we change the orthonormal basis. The 
Slater determinant is precisely this wedge product written in the position basis. There 
is a one-one correspondence between states of the Slater determinant type and points on 
the Grassmannian; there are no longer any spurious variables in the problem. The ground 
state just corresponds to the minimum of E^^(P) on the Grassmannian. 

This point of view can be used perhaps to prove the existence of extrema ( critical 
points) for the Hartree-Fock energy. We make a digresion to point out a connection to 
Morse theory. Gvmi'H) is the well known model for the Classifying Space [21] of U{m). 
So its homotopy type is particularly simple. Also, the cohomology of GrmCH) has a simple 
description; any element can be written as a linear combination of wedge products of the 
forms 

uj2i - tTP{dPf' for i - 1, 2 ■ ■ ■ < m. (179) 

Furthermore there are no relations among these generators. (There would be some relations 
among these generators if were considering instead the Grassmannian of a finite dimen- 
sional vector space.) The odd cohomology groups vanish. The number of generators of 
H^'^ is equal to the number of ways /cwritten as the sum of the numbers 1, 2, ■ ■ -to. Thus 
the generating function of the Betti numbers ('Poincare polynomial') is 

oo m ^ 

p{t) := J2 ^^Hn^ = n jT^j^) (ISO) 

p=0 i=0 ^ ' 

If h and G2 are reasonable operators on Ti and A^Ti, £ will be a differentiable function 
on Gvmi'H). Then Morse theory would guarantee the existence of extrema for £. The 
minimum would be the Hartree-Fock approximation to the ground state; the other extrema 
of finite index will represent excited states. This approach is not possible if we think of the 
HF energy as a function of tta.For, the space of such orthormal frames on Ti (the Stiefel 
manifold of Ti) is contractible; Morse theory would be trivial. There should also be, in 
general, a continuum in the case of atomic physics, corresponding to the scattering states. 
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One can enclose the system in a large box so that all states are discrete and then the 
energy function will only have isolated critical points. 

The weak Morse inequalities say that the number of critical points of index p, Np is 
greater than or equal to dimHP. Thus N2k is greater than or equal to the number of ways 
in which k can be written as a sum of the numbers {1, ■ ■ ■ ,Tn}. One can see that this 
agrees with the counting based on the interpretation of the excited states of an atom in 
terms of electron-hole pairs. The index is the number of directions in which the energy of 
a state decreases, which is related to the number of states with energy less than the given 
state. 

2. Second Quantized approach to Hartree-Fock Theory 

It is clear that mean field theory is a sort of semi-classical approximation to the atomic 
system: the quantum fluctuation are small. So there must be some sense in which Hartree- 
Fock theory describes a classical theory whose quantization gives the atomic hamiltonian. 
However this cannot be the conventional classical theory of particles moving in a Coulomb 
potential, since the atom has no stable ground state in that approach. There must be some 
other classical theory, with a stable ground state, whose quantition also leads to the atomic 
hamiltoinian. We will show that this is a system whose phase space is the Grassmannian, 
and that time-dependent Hartree-Fock theory is equivalent to hamiltonian dynamics on 
this space. 

To begin with let us rewrite the problem in second quantized language. Define the 
creation-annihilation operators associated to the Hilbert space 7i: 

[J{x),a{y)]+ = d{x-y); [a^x),a\y)]+=0; [a{x),a{y)]+ = 0. (181) 

There is a representation of this algebra on the Fock space JF = ©^^QA"^(7i). The vacuum 
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is the state containing no electrons 

a(x)|0>=0. (182) 

The number operator is 

iV^/«t ,,)„,,)., (183) 

and its eigenspace with eigenvalue m is A'^(7Y). The elements of A"^(7Y) are of the form 

\^ >— ip{xij- ■ ■ jXm)ci (xi) ■ ■ ■ a' {xm)dxi- ■ ■ dxm\0 > ■ (184) 

In this language, the hamiltonian is the operator on JF 

H— a^ {x)ha{x)dx + / dxdyG{x — y)a' {x)a^ {y)a{y)a{x). (185) 

The relation to the earlier language is, 

H\tl; >= \H^ > . (186) 

Now suppose we consider a generalization of this atomic physics problem where each 
electron has another quantum number ('color') a = 1, ■ ■ ■ , Nc- We have the operators 
satisfying 

[at-(x), a^{y)]+ = S'^S{x - y); [at«(x), a'ff{y)]+ = 0; K(x), a^(y)]+ = 0. (187) 

There is again a representation of this on the Fock space. But we will allow as physical 
states only those annihilated by 

Q;^ = dx[J'^{x)af3{x) - —a^'^{x)a^{x)]. (188) 

J -^ 'c 

These operators generate an SU{Nc) symmetry. Let us denote this 'color invariant' sub- 
space by JFq: 

\t(;>eJ^o ^^ Q'^\ij>=0. (189) 
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Generalize the hamiltonian to the operator, 



1 ^t. 1 



H =—— / a'°'{x)haci{x)dx -\ — r / dxdyG{x — y) 

Nc J 2N^ J 

[a^'^{x)a^'^{y)af3{y)aa{x) + a^"{x)a^^{y)aa{y)af3{x)]. 

This is a generahzation of the hamiltonian familiar from atomic physics. It is clear that 
if Nc = 1, the color singlet condition on the states becomes trivial and the hamiltonian 
reduces to the previous one. This theory of colored fermions is of physical interest only for 
Nc = 1. Yet, we will show that Hartree-Fock ( mean field) approximation is the large Nc 
limit of this theory. Expansion in powers of -^ will yield the expansions around mean field 
theory. Although A''^ = 1 in the physical case, we know that this -^ expansion is in fact 
a good approximation method in atomic physics. Thus we will see that the conventional 
mean field theory can be understood in terms of a 'replica trick'. 



To see this, introduce the operators 



J_„t. 



P(x.y)^j^al'{x)a,{y). (190) 



' c 



They satisfy the commutation relations 

[P{x, y), Piz, u)] = -^ [6{y - z)P{x, u) - 6{u - x)P{z, y)] (191) 

of an infinite dimensional unitary Lie algebra. * The trace J P{x,x)dx is the number 
operator divieigenvalues of this operator are positive integers; let J-'om be the set of all 
states with eigenvalue m. 

The operators P{x,y) forms a complete system of observables on J^orn- i-e., any 
operator that commutes with P{x, y) for all x,y is a, multiple of the identity. This follows 



* Strictly speaking we should regard the Lie algebra as consisting of the smeared objects 
J dxdyP{x, y)K{x, y) where K is the Kernel of a compact operator. Then there will be no 
divergence problems in the representation theory of this algebra 
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from Schur's lemma if we can see that the representation of the above Unitary algebra on 
J-'om is irreducible. In fact the irreducible representations of the Lie algebra of compact 
selfadjoint operators (and the corresponding group of unitary operators that differ from the 
identity by a compact operator) has been worked out by Kirillov. The representations are 
classified by Young tableaux, just as in the finite dimensional representation theory. Now 
it is straightforward to verify that the representation on J-'om is the one with a rectangular 
Young diagram with Nc columns and m rows. 

Although the P{x^y) form a complete set of observables on Tom, they are not all 
independent of each other; they satisfy a quadratic constraint. To see this, first define the 
normal ordered product of two P's: 

: P(x, z)P{y, u) := P(x, z)P{y, u) - ^5{y - z)P{x, u). (192) 

In terms of the creation-annihilation operators, this means that all the creation operators 
stand to the left of the annihilation operators: 

: P{x, z)P{y, u) := -^J''{x)a^^{y)ac.{z)af3{u). (193) 

c 

The normal ordered product will have finite matrix elements in the limit z ^ y. 

Between a pair of elements |'0 >, |V'' > of J-'om, these operators satisfy the identity 

: P{x,y)P{y,u) : dy\ij >=< ^'| ( 1 - — j P{x,u)\ij > . (194) 

The proof is to note that 

< ^'1 / : P{x,y)P{y,u)dy : \^ = — < ij'\o}'^{x)o}'^{y)aMapiu)\^ > dy 

= -^ / < i^'\a^"{x)af3{u)a^'^{y)ac,{y)\iJ > dy 



Tfl 

1-^) <^'\P{x,u)\i^> 
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We have used 



rfyat^(y)a«(y)|V >= -^^f f dya^^{y)a-,{yM >= mS^\^ > (195) 



for ip E J^Or 



Now we can reformulate the theory entirely in terms of the color singlet variables P. 
The hamiltonian is 

- / dxdyG{x - y)[. P{x,x)P{y,y) : - : P{x,y)P{y,x) :]. 
The normal ordering is necessary to avoid self-energy terms. 

Now we see that with our definitions, Nc appears only as a coefficient of the commuta- 
tion relations of the -P(x, y). Thus -^ plays the role of h; it determines the uncertainty in 
measuring the different components of P{x, y) simultaneously. There are certain states in 
the space J-'om that minimize this uncertainty; they are the analogue of the minimum un- 
certainty wavepackets (coherent states) of the usual canonical quantization. An invariant 
measure of the uncertainty ( or the size of quantum fluctuations) is [3] 

/• . . 

[< i^l : P{x,y)P{y,x) :\ip > - < tlj\P {x , y)\ip X ip\P{y,x)\ip >]dxdy. (196) 

The states that minimize are the highest weight states of some basis. Explicitly, they are 
of the form 

\u>= dxi---dXmUi{xi)---u^{xjn)B^xi)-- -3^x^)10 > (197) 

where 

B'f{x)=J\x)---af^^{x). (198) 

Clearly B' (x) createSince the quantum fluctuation around the expectation values vanish 
in these states, they will have a good classical limit. 
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In the limit Nc = 1 these are precisely the states represented by the Slater determi- 
nant. More generally, they describe an embedding of the Grassmannian GrmCH) into the 
projective space V{J-'om), generalizing the Plucker embedding. We can now calculate the 
expectation value of P: 

<u\P{x,y)\u>=J2<i^)My) ■=Pix,y). (199) 

a 

This is just the projection operator to the subspace spanned by the Ua- Furthermore, the 
expectation value of the hamiltonian is 

<u\H\u>= [dx[( ^ f^:!'^^ l +V{x)P{x,x)] + 
J V dxdy J^^y 

- / dxdyG{x - y)[P{x,x)P{y,y) - P{x,y)P{y,x)]. 

Note that so far we have not made any approximations: these are the exact expectation 
values. Since these Slater-type states are not the exact eigenfunctions of H, time evolution 
will take them into more complicated states. 

However, in the large Nc ( semiclassical) limit the deviations will be small and we will 
be able to decsribe the system completely in terms of the classical vraiable P{x,y). The 
classical phase space is the set of all projection operators of rank m; i.e., the Grassmannian 
Gvmi'H). The commutation relations of the variables P tend to Poisson brackets of the 
classical variables: 

{P(x, y),P{z, u)} = i[5{y - z)P{x, u) - 5{u, x)P{z, y)]. (200) 

These are precisely the ones that follow from the standard symplectic form on the Grass- 
mannian: 

uj = -UrP{dPf (201) 

8 

Thus we see that the limit A'^c ^ oo of our generalized atomic system is a classical 

56 



dynamical system with Gvmi'H) as phase space, u as symplectic form, and 
i?^^(P)=/c^x[(^?^l^)^ +y(x)P(x,x)] + 

- / dxdyG{x - y) [P(x, x)P{y, y) - P(x, y)P{y, x)] 

as the hamiltonian. In particular, the ground state is determined by the minimum of E^^ . 
Also the excited states of the system correspond to stationary points of higher index, as 
anticipated earlier. 

Our derivation shows that time dependent Schrodinger equation tends in the limit 
A^c — ^ oo to the equations 

^^|l^ = {P(x,y),i?^^(P)} (202) 

as the equations of motion. This could be an interesting approach to calculating the 
scattering amplitudes of atoms by other atoms. The small oscillations around a stationary 
point will have some charecterictic frequencies. These will determine the correction to 
order -^ of the energy levels. It is clearly possible to develop a systematic semiclassical 
expansion in -^. 

Thus we have a classical system whose quantization gives the atomic hamiltonian. In 
fact, quantization can be performed by finding an irreducible representation of the algebra 
above. The parameter N^ ( effectively ^) must be a positive integer in order that a highest 
weight representation exist. The physically interesting value happens to be A''^ = 1. The 
success of Hartree-Fock theory shows that even in this extreme case, the -^ expansion 
gives reliable results. 
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